Structure of the flux lines lattice in NbSe2'. equilibrium state and influence of the 

magnetic history. 



A. Pautrat, 1 M. Aburas, 1 Ch. Simon, 1 P. Mathieu, 2 A. Brulet, 3 

C. D. Dewhurst, 4 S. Bhattacharya, 5 and M. J. Higgins 6 

'Laboratoire CRISMAT, UMR 6508 du CNRS, ENSICAEN et 
Universite de Caen, 6 Bd Marechal Juin, F-14050 Caen 4, France. 
2 Laboratoire Pierre Aigrain de I'Ecole Normale Superieure, UMR 8551 du CNRS, 
associee aux universites Paris 6 et 7, 75231 Paris CedexS, France. 
Laboratoire Leon Brillouin, 
CEA Saclay, 91191 Gif/Yvette, France. 
,f Institut Laue Langevin, 6 rue Jules Horowitz, Grenoble, France. 
5 Tata Institute of Fundamental Research, Mumbai-4 00005, India. 
''Department of Physics, Rutgers University, Piscataway, New Jersey 08855, USA. 

We have performed small-angle neutron scattering (SANS) of the flux line lattice (FLL) in a 
Fe doped NbSe2 sample which presents a large peak effect in the critical current. The scattered 
intensity and the width of the Bragg peaks of the equilibrium FLL indicate an ordered structure 
in the peak effect zone. The history dependence in the FLL structure has been studied using field 
cooled and zero field cooled procedures, and each state shows the same intensity of Bragg scattering 
and good orientational order. These results strongly suggest that the peak effect is unrelated to a 
bulk disordering transition, and confirm the role of a heterogeneous distribution of screening current. 
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The peak effect in the critical current of a type II su- 
perconductor is a long standing problem in the field of 
superconductivity. Since its observation [3, 0] , no clear 
consensus has emerged as to its origin. Most models are 
based on the Pippard idea [|| who attributed the peak ef- 
fect to the loss of rigidity of the flux lines lattice (FLL) , 
which falls close to B c2 more rapidly than the pinning 
strength due to inhomogeneities. Larkin and Ovchinikov 
(LO) have treated the complex problem of the elastic 
FLL in the presence of random pinning centers Q . The 
main issues are the statistical calculations of a volume V c 
over which the lattice remains correlated and the corre- 
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sponding pinning force per unit volume F p oc V c . In 
addition to the intense theoretical interest, there is the 
possible link between the FLL bulk disorder and the crit- 
ical current. Other interpretations of the critical current 
deal with dominant pinning of flux lines by the surface 
0, 0], without significant role of FLL bulk order. Re- 
cently, the explanation of the peak effect as a genuine 
phase transition has been addressed both theoretically 
and experimentally, and the peak effect is now considered 
as a classical example of an order-disorder transition in an 
elastic system [7( . A superconductor with a peak effect is 
remarkable by its anomalous transport properties. They 
have been tackled by quite thorough series of measure- 
ments in NbSe2, with a large focus on the metastable 
properties Two coexisting macroscopic FLL states 
with different pinning strengths were observed Q and 
the kinetics between these two states can explain most 
of the anomalous transport properties. Paltiel et al have 
also shown that the critical current in the peak effect zone 
is heterogeneous and more important close to the edge of 



the samples [l(| ■ However, the technique of Hall probe 
does not allow concluding on the structure of the FLL. 
The genuine nature of FLL states with the two different 
critical currents remains unknown, even if one is usually 
assumed to be much more disordered following LO ap- 
proach However, this latter assumption needs to be 
experimentally controlled using a structural probe: in- 
deed, a correlation between the critical current and the 
FLL bulk order can not be a priori assumed, as shown 
by different Small Angle Neutron Scattering (SANS) ex- 
periments 
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1121 Il3j . When no change in FLL order 
is observed whereas the critical current is strongly mod- 
ified [lij ]. an interpretation could be that the pertinent 
structural distortions are at a too small scale to be re- 
solved. However, when important changes in the FLL 
order are seen without any modification of the critical 
[HI, it shows more clearly that the related 
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FLL order is not pertinent for the critical current and 
the dominant pinning mechanism was shown to be at the 
surface. SANS in the presence of a peak effect has been 
mostly used in Nb. Measurements show differences in the 
FLL structure after field cooling and zero-field cooling 
and it was proposed to arise from superheating and su- 
percooling of metastable states as expected for a genuine 
first order (melting) transition [3, EH • Other authors 
proposed that the results would be better understood in- 
volving a crossover into a disordered state [l6j]- The ab- 
sence of intrisinc (thermally driven) melting transition in 
Nb is supported by the persistence of the Abrikosov lat- 
tice up to B C 2 in a purest sample which does not present 
the peak effect [17| . Fe doped NbSe2 systematically dis- 
plays a marked and extended peak effect, and then was 
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studied as a school case using transport measurements 
Q. In addition, no correlation was found between the 
topology of the FLL structure and the enhancement of 
the critical current in NbSe2 using the decoration tech- 
nique (l8j . This makes this sample particularly inter- 
esting for a SANS study of equilibrium and metastable 
FLL, in order to clarify their nature and the link with 
the pinning properties. 

Our sample is a large single crystal of 2H — NbSe2 
doped with 100 ppm of Fe (T c — 5.80K), which was 
used in [3]. SANS experiments were performed at 
ILL (Grenoble-France) using D22 and Dll spectrome- 
ters. The neutron beam was parallel to the magnetic 
field and no external current was applied. In this ge- 
ometry, the higher angular resolution is obtained on the 
rocking curves which measure correlation along the flux 
lines. The neutron wavelength was A — 8 and 10 A with a 
spread AA/A « 10%. Depending on the diffraction vector 
Q of the FLL Bragg planes and in order to optimize the 
intensity, different collimation and sample-detector dis- 
tances were used. The beam was collimated in the middle 
of the sample. Background scattering, obtained from the 
sample in the normal state (B > B C 2) can be subtracted 
to reveal the 2D FLL image (fig.l). All the measure- 
ments were made at a temperature T = 2K, using field 
cooling (FC: magnetic field applied at T > T c ) or zero 
field cooling (ZFC: magnetic field applied at T = 2K) 
procedures. 

A crucial point for a quantitative interpretation of the 
results is the value of the upper critical field B C 2 . It was 
measured using different techniques (magnetization, re- 
sistivity, V(I) curves) on a small crystal cleaved from the 
large one. V(I) curves were also performed on the whole 
crystal used for SANS measurements. All measurements 
show a sharp increase of the critical current after ZFC 
at B — 1AT and a second critical field B C 2 ~ 2.55T at 
T = 2K, as shown in fig.l (ZFC procedure avoids much 
of metastability when a peak effect is present). We ob- 
serve also strong evidences of surface superconductivity 
for B > B C 2- One can estimate the field where the peak 
effect occurs as the midpoint of the abrupt increase of 
I C (B). This gives -Bpeafc ~ 1.5T, i.e . B mn ,k/ 'B r % ~ 0.6 
in agreement with previous reports (l9| . First series of 
SANS measurements relate to the structure of FLL in 
the area of the peak effect, to see if there exists a transi- 
tion toward a disordered phase. Such measurements were 
made for applied magnetic fields from B = OAT to 2T, 
using the ZFC procedure to directly compare with the 
critical current measurements shown in fig.l. For fields 
larger than B > 2T, the scattered intensity from the 
FLL became too small to measure with counting times 
of 4 Hours on D22, but this value is however much higher 
than B pea k- Such measurements always show a hexago- 
nal lattice with sharp Bragg peaks and a diffraction vec- 
tor Q corresponding to the applied magnetic field. For 
each field value, we have measured the full rocking curve 



to get the scattered intensity of the (1,0) Bragg peak of 
the FLL. Each rocking curve was fitted by a Lorentzian 
(fig. 2) in order to extract the integrated scattered inten- 
sity and rocking curve width. 

The intensity for a (hk) reflection is 



Jwt - 2TrVcf > ( 1 /4) 2 X 2 n /cj ) o | F hk | 2 /Q hk 
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where V is the illuminated sample volume, cj> the incident 
flux, 7 the magnetic moment of the neutron, and <po the 
flux quantum. For a quantitative analysis of the scattered 
intensity at high fields, the Abrikosov limit can be used 
in the range 0A.B C 2 < B < B C 2 20], where the Fourier 
factor of the hexagonal vortex lattice can be expressed 
as: 
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where v = (h 2 + k 2 + hk) and n Q .M = (B - 
B C 2) I f3a(2K 2 — 1) according to the Abrikosov limit ((3 a — 
1.16 for the hexagonal lattice, re is the Ginzburg-Landau 
parameter) . 

After normalization to the incident neutrons flux, F\q 
can be extracted from equation [T] We choose to fix 
B c2 = 2.55T which was robustly determined, and only 
one parameter re is adjustable. i*io compares favorably 
at high fields with the Abrikosov form of the intensity 
with re 12 (fig. 2), giving A = re.£ = 136nm with 
£ 2 = 4>o/{2ttB C 2) « 128nm 2 . To describe the data at 
the lowest fields, we used numerical corrections from 
Ginzburg-Landau solution as described in 0, [5l| . Since 
it converges to the Abrikosov limit for B/B C 2 > 0.4, 
there is no supplementary fitting and the agreement with 
the experimental data is satisfactory (fig. 2). The val- 
ues of A and £ are consistent with the literature, con- 
sidering that previous estimates of A are quite scattered 
over the broad range 70 — 250nm ( (22j | and references 
herein). The full- width-half-maximum (FWHM) of the 
rocking curves was analyzed with the correction of in- 
strumental resolution, mainly due to wavelength spread 
and the beam divergence. It can be analytically esti- 
mated as a first approach using a Gaussian approxima- 
tion [23|]. For a peculiar set-up, we find for example 
FHWM (resolution) = OAGdeg for B = OAT, to com- 
pare with FW H M (rocking — curve) = 0.22deg, and an 
intrinsic with Aoj can be extracted [24j |. Finally, a lon- 
gitudinal correlation length d ss ao/irAuj is deduced. As 
shown in figure 3, we have found that d ~ 11 ± Aum is 
roughly constant for all field values (the crystal is 200 um 
thick). Note that, since FWHM is close to the instru- 
mental resolution, the value of d has rather the meaning 
of a lower bound 25j. This value is however already far 



larger than the one deduced from a criteria used to in- 
dicate the appearance of dislocations in the FFL (the 
order state is expected to be unstable for d < 20ao (2r3 | 
and we measure d > 200ag). Due to the poor resolution 
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in the positional correlation function in our SANS geom- 
etry, we can not distinguish between a perfect lattice of 
flux lines, a Bragg glass state with algebraically decay- 
ing translational order 27| , or a multi-domains structure 
which fractures at large scale [28j]. Detailed investiga- 
tions of correlation functions of the FLL using a high 
resolution geometry indicate some large scale fracturing 
of FLL even in best quality Nb samples [29| • We summa- 
rize the most important result of this part: our analysis 
of FLL Bragg peaks stabilized after ZFC shows that the 
FLL order is unchanged and robust far inside the peak 
effect where the critical current has strongly increased. 

In NbSe2, FC states show different properties than 
ZFC states [8(. In particular, they present a large and 
metastable critical current and are macroscopically het- 
erogeneous, even for magnetic fields smaller than B pea k 
[1, EDI ■ This is a quite unusual case, since FC leads usu- 
ally to a more perfect lattice in samples without the peak 
effect. Comparison between FC and ZFC scattering in- 
tensity was here tricky. Rocking the sample through the 
expected Bragg angle after FC gives a scattered intensity 
significantly reduced. One might conclude that much of 
the intensity was incoherently scattered from a strongly 
disordered supercooled state. However, as noted in [301 ] 
when measuring CeRu2, if the Bragg peaks are strongly 
broadened, all the intensity is not recovered for restricted 
angular scans. This is also the case if the Bragg angle has 
changed. Measurements performed with large enough 
angular scans show that most of the intensity in the FC 
rocking curve was surprisingly outside the regular Bragg 
conditions (fig. 4). The shape of the rocking curve is dif- 
ficult to extract because it emerges only slightly from 
the background. However, an analysis with two simi- 
lar Bragg peaks rotated from the theoretical Bragg angle 
is most satisfactory. Even if two peaks in the rocking 
curve indicates a very peculiar FLL structure, the sum 
of integrated intensity of the two peaks corresponds to 
the intensity expected for the regular Abrikosov lattice 
measured after ZFC (fig. 2) [3l[. The Q value is also un- 
changed and its radial width is resolution limited. Since 
the intensity in the RC is preserved, all the "disorder" 
should be along the longitudinal direction, without af- 
fecting the orientational quality of the FLL. This is con- 
sistent with measurements of the orientational width of 
the Bragg peaks in the plane of the detector. They re- 
main resolution limited (Aip = 13Hdeg in the azimuthal 
direction) for the FC case (fig. 5). Due to this relatively 
poor azimuthal resolution, we do not know if the FC FLL 
presents an orientational order as good as the ZFC FLL. 
But clearly, the FC FLL is not a supercooled glassy-like 
state which should lead to a degenerated orientational 
order (a ring of scatterin g) 11411 . The double peak struc- 
ture was first observed in [13( . It implies that the FLL is 
tilted from the magnetic field direction due to magnetic 
field components in the plane. A non symmetric dis- 
tribution of superficial screening currents was proposed 



as a qualitative interpretation [13j. The new informa- 
tions obtained here (the preserved intensity of the whole 
FLL after FC, and the absence of bulk transition in the 
ZFC FLL when crossing the peak) are consistent with 
this scenario. We have tried the usual "shaking" pro- 
cedure consisting in a variation of the magnetic field in 
the form of a sinusoidally damped oscillation, in order to 
improve the FLL spatial order [32|]. The effect was only 
marginal, likely due to the non-trivial nature of the FC 
FLL structure. Contrarily, a thin Bragg peak centered 
at the theoretical Bragg angle is recovered when pass- 
ing a large dc overcritical current in the sample [l3j. A 
possible reason of the negligible effect of the "shaking" 
is that the screening currents are too large to be affected 
by the currents induced by the field oscillation. Another 
one is that the screening currents have to be polarized 
in one direction to recover an homogeneous distribution 
and a regular FLL. Note that our results appear consis- 
tent with the edge contamination model [Io| and with the 
observation of two coexisting macroscopic vortex states 
[9(, to the extent that the supercooled state with the 
high critical current stabilizes surface currents only and 
is unrelated to a bulk transition. Such non homogeneous 
distribution of surface currents is also consistent with the 
complicated current flow necessary to induce the inter- 
face dynamics between the coexisting FLL states [9( . The 
origin of two competing surface currents is unclear, but 
may arise from a coexistence of two separate mechanisms 
contributing to the surface transport current in the mixed 
state. Such a scenario was proposed some years ago from 
angular magneto-transport measurements coupled with 
surface treatments in conventional type II superconduc- 
tors (33| , and deserves a peculiar attention. 

In conclusion, we have studied the FLL structure in a 
Fe doped NbSt2 crystal, which shows a large peak effect 
in the critical current. The ZFC FLL is always ordered. 
The FC FLL is notably different, but does not present 
the characteristics of a bulk, glassy-like, disordered state. 
These results are extremely difficult to conciliate with a 
fracturing or melting transition of the FLL usually pro- 
posed as the cause of the peak effect and calls for al- 
ternative interpretations. A direct consequence of our 
SANS study of FLL structures is that the peak effect may 
correspond to an anomalously large surface or boundary 
current developing in the peak. It can be frozen over 
macroscopic zones of the sample and leads to the two 
states mixture responsible for the metastable transport 
measurements. 
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FIG. 1: Left: 2D pattern of the ZFC FLL (B = OAT, T = 2K) 
obtained by rocking the sample through the Bragg condi- 
tions and after a subtraction of the normal state background. 
Right: critical current Ic(B) measured after ZFC at T = 2K 
in superfluid Helium and showing the sharp increase of critical 
current: the peak affect at B pea k ~ 1.5T. 
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FIG. 2: The full points are the experimental Fourier compo- 
nent Fiq deduced from the scattered intensity Jio of the ZFC 
FLL and compared to the Abrikosov and to the Ginzburg- 
Landau models (B C 2 = 2.5T, A = 120nm). The empty points 
are experimental points corresponding to the FC FLL. In 
the inset is shown the square of the Fourier component to 
illustrate the expected cross-over at B w QA.B C 2 = IT from 
first order Ginzburg-Landau correction to the London model 
at low field (shown as a straight line in this scale) to the 
Abrikosov line. 
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FIG. 3: Correlation length d deduced from the rocking curve 
width of a first order Bragg peak as function of the magnetic 
field (T = 2K) after ZFC. In the inset is shown a typical 
rocking curve better fitted with a Lorentzian (B — 1.25T). 
No background is subtracted. 
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FIG. 4: Top: rocking curve measured after a zero field cooling 
for B = OAT. The fit is a Lorentzian with FWHM = 0.13 ± 
O.Oldeg. Bottom: rocking curve measured after a field cooling 
for B = OAT with the same scale as the top figure. The 
inset is a zoom showing two Bragg peaks. The fit is with 
two Lorentzians of FWHMi = 0.7 ± 0.2deg and FWHM 2 = 
1.0 ± O.ldeg which are roughly symmetric with respect to the 
theoretical Bragg Angle. No background is subtracted. The 
Q value is the same for FC and ZFC FLL. 
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FIG. 5: Azimuthal scan for FC and ZFC states (B = OAT). 
Despite ui rocking curves are very different, the orientational 
order remains resolution limited as demonstrated by the mean 
value Atp of the peaks FWHM (ZFC:At/j = 12.98 ± 0.50deg, 
FC:14.09 ± 1.38deg). For the FC case, all the 6 peaks are 
observed without performing full scans for all the (h,k) reflec- 
tions because the lattice is in Bragg condition for an extended 
angular range. 



